The static and dynamic conductivity of warm dense Aluminum and Gold calculated 

within a density functional approach. 
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The static resistivity of dense Al and Au plsmas are calculated where all the needed inputs are 
obtained from density functional theory (DFT). This is used as input for a study of the dynamic 
conductivity. These calculations involve a self-consistent determination of (i) the equation of state 
(EOS) and the ionization balance, (ii) evaluation of the ion-ion, and ion-electron pair-distribution 
functions, (iii) Determination of the scattering amplitudes, and finally the conductivity. We present 
data for for the static resistivity of Al for compressions 0.1-2.0, and in the the temperature range 
T — 0.1 - 10 eV. Results for Au in the same temperature range and for compressions 0.1-1.0 is 
also given. In determining the dynamic conductivity for a range of frequencies consistent with 
standard laser probes, a knowledge of the electronic eigenstates and occupancies of Al- or Au 
plasma become necessary. They are calculated using a neutral-pseudoatom model. We examine a 
number of first-principles approaches to the optical conductivity, including many-body perturbation 
theory, molecular-dynamics evaluations, and simplified time-dependent DFT. The modification to 
the Drude conductivity that arises from the presence of shallow bound states in typical Al-plasmas is 
examined and numerical results are given at the level of the Fermi Golden rule and an approximate 
form of time-dependent DFT. 



PACS numbers: PACS Numbers: 52.25.Mq, 05.30.Fk, 71.45.Gm 
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I. INTRODUCTION 



The static conductivity of matter in the plasma state 
can be calculated with some confidence, at least for 
a number of "simple" plasmas, using an entirely first- 
principles approach p], 0- Comparisons with experi- 
ments are available for a wide range of conditions, at least 
for Al plasmas |3| ■ Another transport property of great 
interest is the frequency-dependent conductivity o~(lj). 
Optical probes provide a convenient diagnostic tool for 
this plasma conductivity, and common laser-probe wave- 
lengths, e.g., 300-800 nm, become the window of interest. 
In fact, many experiments on laser-plasma interactions 
related to inertial fusion have concentrated on frequency 
tripled (with the fundamental u> at 1/im from glass lasers) 
light, 3w, 351nm, while there is also current interest in 
the 2u> (527 nm) regime as an option for indirect-drive 
ignitionQ. Normally, if the probe frequency is less than 
the plasma frequency, the probe photons fail to enter 
the material and only weakly ionized systems can be ac- 
cessed. However, the recent development of extremely 
thin (" nanoscale" ) plasma-slab techniques 5] , known as 
idealized-slab plasmas |6j has made it possible to study 
dense plasmas with low-frqeuency probes, both in trans- 
mission and in reflectivity. 

Density-functional calculations are of two types. The 
first type depends heavily on quantum MD (QMD) sim- 
ulations, e.g, Car and Parrincllo treat only the electrons 
via Kohn-Sham theory, while the ion subsystem, explic- 
itly represented by a convenient number N (of the or- 
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der of ~ 32-256) of ions, is made to evolve in time and 
an average over millions of configurations is taken. In 
full Quantum Monte Carlo (QMC) simulations [7J,|8j, even 
Kohn-Sham theory is not used and hence full QMC is not 
practical for typical plasma problems. The second type 
of DFT is typified by our approach where both electrons 
and ions are treated by DFT, so that both subsystems are 
described by two coupled "single-electron" and "single- 
ion", Hartree-like Kohn-Sham equations. The many- 
body effects (i.e, many-electron and many-ion effects) are 
included through the exchange-correlation or correlation 
potentials. This allows an enormous simplification in the 
numerical work involved. These equations may be fur- 
ther reduced to extremely simpilified Thomas-Fermi ap- 
proaches yielding approximate results for plasma proper- 
ties which are within an order of magnitude of the more 
refined results, even in unfavourable cases 9J. 

The general optical conductivity problem is essentially 
the same as that of opacity calculations for warm dense 
matter^cj. However, here we approach it from the static- 
conductivity (u> — 0) regime, and take account of free- free 
processes as well as shallow bound states which may ex- 
ist. The static-conductivity provides an evaluation of the 
dynamic conductivity near u> = via the basic Drude for- 
mula which assumes a constant relaxation time r, taken 
to be the static collision time r(0). Hence our task can 
be stated as follows: 

1. Calculation of the static resistivity. This involves 
the following steps: 

(a) Calculate the Kohn-Sham atom in an electron 
gas of density n. 

(b) Use the Kohn-Sham results to form pseudopo- 
tentials Vi e and scattering cross sections at the 
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given density p and temperature T. 

(c) Use the Vi e to form pair potentials and pair- 
distribution functions. Here the Kohn-Sham 
equation for the ion-subsystem can be ap- 
proximated by some form of the hyper-netted- 
chain equation (HNC) with bridge terms. 

(d) Calculate the EOS, ionization balance etc., to 
obtain effective ionic charges Z , electron den- 
sity n = pj Zetc, and self-consistently, repeat- 
ing from item (a). 

(e) Calculate the static resistivity and the static 
relaxation time r(0) using, e.g, using a Ziman- 
type formula valid for strong coupling and fi- 
nite T. 

2. Use the energy-level structure of the Kohn-Sham 
atom for the p and T regimes of interest to set 
up bound-bound, or bound-free processes which fall 
within the range of frequency uj considered. Correc- 
tions may be necessary since the Kohn-Sham the- 
ory does not provide good excitation energies. 

3. Construct the dynamic conductivity a°(uj). This 
yields results equivalent to the Fermi golden rule. 

4. Extend the calculation of cr(uj) u sing time- 
dependent density-functional theory flfj .Include 
the coupling of electronic transitions to ion- 
dynamics. 

While some aspects of this program can be carried 
out, it is fair to say that a complete, consistent theory 
of the transport properties of such many-body systems 
as posed by warm dense matter, or indeed, even the 
static properties as embodied in the equation of state 
for some regimes of density and temperature are open to 
debate, even for well studied systems like aluminum [11J, 
and hydrogen^^. The objective of the present study is 
to calculate the dynamic conductivity of Al and Au plas- 
mas for a number of compressions and densities where 
the Drude theory may need correction due to the pres- 
ence of shallow bound states. Such boundstates ionize 
when the compression is changed, and produce distinc- 
tive changes in transport properties. Sharply rising static 
resistivities under a change of compression is a common 
feature of some of the theoretical calculations shown here. 
However, the establishment of a genuine phase transition 
requires more careQ. The cr(w) of expanded liquid met- 
als and plasmas show^l effects arising from clustering 
and excitonic effects, as the metal-insulator transition is 
approached. These excitonic effects are not important in 
dense systems. While Al has been an object of exten- 
sive study, recent experiments in the warm dense matter 
regime have focused on gold targets 0, 0|. Here we 
present numerical results for Al and Au for several com- 
pressions and temperatures. 



II. STATIC RESISTIVITY 

We use atomic units (Hartree=l a.u., the Bohr radius 
ao = 1, with |e| = h = m e = 1). The atomic unit of 
resistivity, given by hao/e 2 has a value of 21.74 /iSlcm. If 
the equilibriation of the electron distribution perturbed 
by the applied electric field is governed by a relaxation 
time r, the conductivity a is given as: 

A mean free path l m f p — < v > /r, where < v > is some 
characteristic mean velocity, is often introduced. If elec- 
trons were classical point particles, then it is evident that 
Impf cannot be smaller than the mean separation between 
collision centers. This is sometimes called the Joffe- 
Regel-Mott rule, and holds well in many semiconductors. 
However, electrons are quantum particles (wave packets) 
with an extension of the order of the thermal de Broglie 
length. Further, r depends on the electron momentum 
k. Thus there are examples where l m pf, obtained by 
some averaging process, is in fact smaller than some esti- 
mated "mean-ion separation" . Although the concept of 
the mean free path is implicit in the Boltzmann-equation 
approach, this is not necessary in the Dyson equation 
which replaces the Boltzmann equation in the quantum 
case. If the mean free path is large, the "particle picture" 
of the electron applies, while if this is comparable to the 
lattice parameter, then we are in the diffraction limit and 
the wave picture must be used. The Dyson equation, 
valid at both extremes, describes the one-particle propa- 
gator which is closely related to the distribution function 
appearing in the Boltzmann equation. The derivation 
of a transport coefficient from the Dyson equation leads 
us to the current-current correlation function. This is 
closely related to various two-body distributions and col- 
lision kernels found even in classical kinetic models. Un- 
fortunately, although formal expressions can be written 
down, their evaluation using Green's functions or related 
methods becomes impractical, especially if bound states 
are present. In effect, Green's-function methods can be 
pushed to, say, second order in the screened interactions. 
Such an approach is sufficient if there are no boundstates 
associated with the scattering potential. Attempts to go 
further rapidly become intractable and useless. 

Our approach is to use a variety of techniques and 
replace the ion-electron interactions by suitably con- 
structed pseudopotentials, or use scattering cross sections 
calculated from phase shifts. This requires a fairly so- 
phisticated non-perturbative description of the ion-ion 
and ion-electron correlations in the plasma. 

A. Description of the plasma using the Kohn-Sham 
equations 

We begin with the bare nuclei and construct the elec- 
tronic and ionic structure of the plasma. The interaction 
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of the nucleus with the electron fluid is a highly nonlin- 
ear process and attempting to treat it using perturba- 
tion theory is unfruitful. Hence we use the Kohn-Sham 
technique, and construct the non-linear charge density 
around the nucleus. The nucleus, together with its charge 
cloud of bound states and continuum-electron states con- 
stitute a neutral object. This neutral object is called 
the neutral pseudo-atom (NPA), following the usage of 
e.g., Ziman and Dagens Thus an important result of 
the Kohn-Sham procedure is the charge-density "pilcup" , 
n(r), around the nucleus that essentially screens the nu- 
cleus. A part of this arises from the free electrons and 
is denoted by rif(r). This n(r) and n/(r) depend on 
the mean density n of the electron fluid, temperature T, 
and the nuclear charge Z . The Kohn-Sham procedure 
leading to the NPA provides the phase shifts suffered 
by the continuum electrons when they scatter from the 
nuclei. These are used for constructing the scattering 
cross sections (or pseudopotentials) which describe the 
electron-ion interaction. The pseudopotential has an ef- 
fective ionic charge Z and it behaves as —Z/r for large r. 
The rapid oscillations of the potential near the nucleus 
is replaced by a weak, smooth core region as the "va- 
lence" electrons do not really penetrate the ion-core. The 
pseudpotentials used in many of the solid-state or molec- 
ular code packages ^3 have the necessary transferability 
and could be quite useful. However, they assume that 
the pseudopotentials would be used within a Schrodinger 
or Kohn-Sham type procedure rather than in a linear- 
response scheme, and hence they cannot be directly used 
within a Ziman-type formula. 

Sometimes, instead of using the all-electron Kohn- 
Sham equation or using a suitably constructed pseu- 
dopotential, the electron-ion interaction is replaced by 
a Yukawa-type interaction (effectively, a Debye-screened 
interaction). The electronic structure is calculated for 
such a static-screened nucleus. This procedure is not 
justifiable since the energies of bound-state electrons 
correspond to very high frequencies at which there is 
no screening. Further, the orthogonality of the contin- 
uum eigenstates and the bound states ensures that there 
is very little penetration of the free electrons into the 
bound-electron region. Consequently, there is very little 
screening of the inner bound states, where as the Yukawa 
potential screens the inner bound states as well. Thus 
calculations using a Debye-like potential in warm dense 
matter is likely to be incorrect irrespective of whether 
the Born approximation, or a T-matrix, etc., were used. 

In the case of Al and also Au, it is possible to construct, 
in many situations, a soft pseudopotential Vi e (q) which is 
weak in the sense that it is possible to recreate the non- 
linear electron-density "pileup" H/(r), obtained via the 
Kohn-Sham equation, by within linear response theory. 
That is, we define the Vi e (q) such that 

n/(r) = -V ie (q) X (q) (2) 

Here x(l) is the electron linear-response function. Un- 
like the transferable pseuodpotentials used in it ab initio 



packages, this Vi e (q) is specific to the chosen Z, n, T, 
and the atomic number Z. It is often convenient to write 
the pseudopotential in the form 

V ie {q) = ZM q v(q), v q = 4ir/q 2 (3) 

where v q is the bare Coulomb potential and M q is a form 
factor. Only a local pseudopotential is used, and this is 
quite adequate for an analysis of the experimental data 
currently available. Since the pseudopotential is weak by 
construction, the ion-ion pair potential can be taken to 
be 

U ll (q)=^Z 2 /q 2 + \V te (q)\ 2 x(q) (4) 

Given Uu(q), the ion- ion distribution function gu{r) 
and the structure factor S{k) can be calculated using 
the Hyper-netted-chain (HNC) equation or its extension 
where a bridge term is included. We note (see below) 
that the HNC equation (or its extension) is in fact the 
Kohn-Sham equation for classical particles ( here, Al 2+ 
or Au 2+ ions ) for certain choices of the ion-ion corre- 
lation potential (there is no exchange potential because 
the ions are classical particles). 

To summarize, by using the Kohn-Sham procedures, 
we have thus obtained the ion-electron pseudopotential 
Vi e (q), the structure factor Su(q), and the charge density 
riie(r), and an effective ionic charge Z which enters into 
the pseudopotential. The phase shifts have been used 
to construct a scattering cross section pj which may be 
used instead of the pseudopotential. Hence we have a 
completely self-consistent procedure for obtaining all the 
relevant quantities starting from the nuclear charge of the 
element. The numerical codes for carrying out these pro- 
cedures are available via the internet, to any interested 
researcher Ha. 



B. Is this a "one-center" approach ? 

To answer this question, we consider the density- 
functional theory of a two component system consisting 
of electrons (density profile n(r) ), and ions, with a den- 
sity profile p(r), with respect to an ion positioned at the 
origin |l9|. Then the Hohenberg-Kohn-Mermin theorem 
states that the free energy F[n(r), p{r)] is a functional of 
the density distributions such that: 

dF[n{r).n{r)]/dn{r) = (5) 
dF[n(r),p(r)}/dp(r) = 0. (6) 

The first of these equations leads to an effective single- 
electron equation, viz., the Kohn-Sham equation where 
the effective potential contains an exchange-correlation 
potential which takes account of many-electron effects. 
The second equation also leads to a Kohn-Sham equation 
which is a classical equation for a single ion. This also 
contains an ion-ion correlation potential which brings in 
the effects of the multi-centered system. It can be seen 
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that this classical Kohn-Sham equation reduces to the 
HNC equation for a certain choice of the ion-correlation 
potential. Further, for "simple" metallic plasmas like 
Al where the pseuopotential is weak, this scheme re- 
lates closely to pair-potential based liquid-metal theory. 
Such a simplification does not hold, for e.g., for hydro- 
gen plasmas|l2j. Since ion-ion, ion-electron and electron- 
electron correlations are included in the theory, it is not 
a single-ion model of the plasma. It is firmly rooted in 
a many-electron, many-ion DFT approach which does 
not invoke the Born-Oppcnhcimcr approximation. Since 
the theory is explicitly based on distribution functions, 
it is manifestly non-local and can be easily implemented 
to be free of electron self-interaction errors. Our ap- 
proach may be contrasted with the Car-Parrinello (CP) 
approach where DFT is used only for the electrons, while 
the ions are explicitly and individually treated by classi- 
cal molecular dynamics. CP avoids the need for an ion- 
correlation potential, but demands a much larger compu- 
tational effort. Some of these efforts, based on MD meth- 
ods (e.g, Ref. 0, 0| ) have confirmed results obtained by 
the numerically simpler methods that we have used. 

C. Extended Ziman formula for strongly-coupled 
electrons and ions. 

The Ziman formula is an application of the Boltzmann 
equation to liquid metals. It was extended to finite tem- 
peratures by a number of authors. |22j| The crux of the 
problem is the evaluation of the collision rate. Com- 
pared to some methods well known in plasma theory (e.g, 
Lenard-Balescu) the collision rate is easily evaluated us- 
ing the "Fermi golden rule" . In this section we briefly 
recapitulate Ziman theory in the language of the Fermi 
golden tulc. In the relaxation-time approximation we as- 
sume that the perturbed Fermi distribution f(k) for elec- 
trons with momentum k relaxes towards the equilibrium 
distribution fo(k) according to the equation: 

_ f(k) - fo(k) 
dt lc ° l r(fc) 1 ' 

Considering an electron scattered elastically from state 
k to state k' , with = \k'\, e k = e k ', the net scattering 
rate is the difference of the two processes (k — > k') — (k' — > 
k). The initial and final densities of states for the k — > k' 
process is f(k) and [1 — f(k')]. Hence the Fermi Golden 
rule gives 

R(k^k') = (2ir/h)Y,\Tkk>\ 2 5(e k - e k ,)f(k)(l- /(*') 
R(k' — > k) = permutation of k with k' etc., 

Since e k = ey , this involves only an angular integration 
and \k\ = \k'\. Since the energy is not changed, the static 
resistivity arises purely from momentum randomization. 
The change in momentum is q 2 — 2fc(l — cos(9)). Here 9 
is the angle between k' and k. This (1 — cos(6)) term does 



not appear in the usual relaxation time which is the time 
between scattering events. Using these rates in Eq.[7| we 
obtain a result for the inverse of the relaxation time. 

-L = 2tt£ 6(e k - e' k )\T kk ,\ 2 (l - cos(9)) (8) 

Here the sum merely indicates the integration over 9. The 
"T-matrix" appearing here describes the scattering of an 
electron by the whole ion-distribution (i.e, not just one 
ion). Given a set of ions at instantaneous positions Ri, 
then the interaction of an electron at r with the whole 
distribution is of the form: 

V(r) = Y J V le (r-R l ) (9) 

The matrix element between the initial state k and the 
final state k', with q = k' — k is: 

V{q) = Y J ZM q v q e 1 ^'. (10) 

Note that 

P q = ^2exp(iq- Rj) (11) 

< PqP - q > = NiS u (q). (12) 

Thus the ion-ion structure factor S(q) and the single-ion 
scattering cross section (or the pseudopotential) from a 
single ion combine to give the full scattering T-matrix. 
The dependence on the structure factor becomes negligi- 
ble for T > 10 eV. The individual scattering cross section 
can be replaced by a single-center T-matrix (to be de- 
noted by ifcfc/ ) obtained from the phase shifts of the NPA 
calculation for a single nucleus. The S(q) is also obtained 
from the pair-potential constructed from the same NPA 
calculation. In Ref. |lj we showed how to avoid the cal- 
culation of the S(q) by directly computing the scattering 
cross section from the whole ion distribution. That is, 
T kk > s not factored into single-center t kk i and the asso- 
ciated structure factor. Such an approach is needed for 
strongly interacting systems where such a factorization 
may not be valid. In this context we note that the resis- 
tivities calculated by us using the "single-center" model 
(ie., t kk >S(q)), and the full ion-distribution model ( T kk > 
for strong-scattering) for H-plasmas were independently 
confirmed by the quantum Monte Carlo simulations of 
Kwon et al(3- However, if the pseudopotential is weak, 
the S(k) and the single-center t{k) may be used. 

Given the inverse relaxation time l/r(fc) for an elec- 
tron of momentum k, or equivalently, l/r(e) for the en- 
ergy, e = k 2 /2, we need to average this over all electron 
energies to obtain a resistivity or a conductivity. The 
averaging used in the Ziman formula leads to a resistiv- 
ity, while a direct application of the Boltzmann equation 
would lead to a conductivity. Thus, 

a = ^<T(e)>, R=^<l/r(e)> (13) 
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The Boltzmann equation shows that the averaging rel- 
evant to the conductivity calculation is such that 

On the other hand, the averaging over the 1/t(c) used in 
the extended Ziman formula for the resistivity is some- 
what different. 

< l/ T >=-J™deg(e)^T(e)- 1 (15) 

Here g(e) is a density-of states factor which is unity for 
free non-interacting electrons. It should be constructed 
from the phase shifts in a strongly scattering environ- 
ment. 

The initial assumption, Eq. [7| was that the modified 
distribution was defined via a relaxation time. A more 
complete approach is to represent the modified part 5f(k) 
as a series expansion in a set of suitably constructed or- 
thogonal polynomials, and obtain a variational solution. 
The set of polynomials appropriate for degenerate (and 
partially degenerate) electrons has been discussed by 
Allen. (23. In the classical limit, such polynomials are the 
well known Sonine polynomials. The usual relaxation- 
type approach is equivalent to a single-polynomial solu- 
tion. This is adequate for dense Al-plasmas, and for the 
range < T < 10 eV studied here. However, this is 
probably not so for Al at 1 /4 of the normal density, or 
for lower densities. We have less experience with Au- 
plasmas to assess the quality of the resistivities for Au 
obtaind here. 



D. Numerical results in the static limit. 

The static resistivities for Al calculated using the above 
methods (and some extensions of it) have been compared 
with experiment by Benage et al.0 Given the uncertain- 
ties in the experiment and the various approximations in 
the theory, the agreement is quite good. However, there 
are a number of difficulties in the calculation. If we con- 
sider an Al-plasma at a compression k=0.25, at T > 2.5 
eV the Al-ion has the bound shells Is, 2s, 2p, and 3s. 
The 3s level becomes increasingly shallow as the temper- 
ature is reduced. Below approximately T ~ 2.5 eV, the 
3s level begins to "evaporate" and becomes free, i.e, the 
3s bound electrons ionize. Even at 2.5 eV, the occupa- 
tion of the 3s level is ~ 0.744. The ionization state Z 
increases from ~ 1.667 to 3. The model where we have 
a single average Z is clearly not satisfactory in such a 
region. The plasma contains an equilibrium mixture of 
several ionization states Z^ with concentrations x\ such 
that 

Z = T,iXiZi (16) 

For such situations (and indeed in general), the concen- 
trations Xi have to be determined from the minimum of 



TABLE I: Resistivity (in /xohmcm) of Al plsamsa as a func- 
tion of the compression k and the temperature in eV, calcu- 
lated within the approach described in the text. 



T| K — ► 


0.1 


0.25 


0.5 


1.0 


2.0 


4.0 


0.10 


1620 


801 


77 


28 


38 


48 


0.25 


1582 


812 


84.1 


30.5 


37.4 


47.0 


0.50 


1527 


847 


92.1 


34.4 


35.8 


45.6 


0.75 


1447 


859 


99.7 


36.7 


34.5 


45.0 


1.00 


1340 


846 


107 


38.7 


34.5 


44.8 


1.75 


1063 


707 


122 


44.4 


34.2 


44.6 


2.50 


844 


597 


135 


50.2 


35.4 


44.7 


3.75 


1196 


448 


148 


60.4 


39.5 


46.0 


5.00 


916 


434 


160 


70.8 


45.1 


48.3 


7.50 


646 


380 


180 


90.6 


57.7 


54.7 


10.0 


544 


379 


194 


108. 


70.1 


61.8 



the total free energy. If n species of ions are possible, 
then n different Kohn-Sham equations have to be solved, 
and n(n + l)/2 ion-ion distribution functions have to be 
determined, and the total free energy has to be calculated 
as a function of n and Xi . The minimum property of the 
free energy yields the equilibrium plasma conditions from 
which the resistivity is calculated, using the n scattering 
cross sections and the n(n + l)/2 structure factors. We 
reported such a calculation for Al in Ref. Q. Our ex- 
perience is that, even with a multi-ion fluid model, the 
self-consistent equations may fail to converge since the 3s 
level (or the 2p level) affects the iterative procedure. Ma- 
jumdar and Kohn have shown that physical properties of 
a system should be continuous across the region where an 
electronic state moves from being a bound state to a con- 
tinuum state[24j. Hence we may calculate the resistivity 
in two adjacent regions separated by a non-convergent re- 
gion, and smoothly join the calculated resistivity across 
the "difficult" region. In Table [H] we show the Kohn- 
Sham eigenvalues of the bound-electron states for Al ions 
in various plasmas. Although the Kohn-Sham eigenval- 
ues do not exactly correspond to the excitation energies, 
they provide an initial estimate which can be improved 
using the methods of time-dependent density functional 
theory, or using self-energy calculations!^. As the 
density decreases and the temperature increases, we in- 
evitably pass through regions of T, k where the problem 
of bound states which hover near ionization becomes im- 
portant. We return to this question in discussing the 
dynamic conductivity of Al-plasma at k = 0.25, and 0.1. 
Table I shows that, for k = 0.1 the resistivity essentially 
decreases with T, while for k = 0.25 the resistivity grad- 
ually increases in value, goes through a plateau-like re- 
gion, and then begins to decrease with temperature. The 
same behaviour holds for the other higher compressions, 
but the table given here does not go high enough in tem- 
perature (for the higher densities) to show the plateau 
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FIG. 1: Resistivity of warm dense Aluminum as a function of 
the compression for several temperatures 



FIG. 2: Resistivity of warm dense gold as a function of the 
compression for several temperatures 



TABLE II: Kohn-Sham energy-level structure for several Al- 
plasmas within the NPA averag-configuration model 



Level 


K = 4 


K = 2 


K= 1 


k = 0.5 


k = 0.25 


K = 0.1 




T=10 












Z 


3.043 


3.0166 


3.0164 


3.0194 


2.6060 


2.2427 


2s 


-5.7016 


-6.3615 


-6.9213 


-7.4089 


-7.9601 


-8.5825 


2p 


-2.9311 


-3.5922 


-4.1529 


-4.6411 


-5.1928 


-5.8159 


3s 










-0.1871 


-0.6168 


3p 












-0.1643 




T=2.5 












Z 


3.0427 


3.0051 


3.0003 


3.0000 


1.6643 


1.6240 


2s 


-5.5247 


-6.1819 


-6.6655 


-7.0257 


-7.4435 


-7.6634 


2p 


-2.7838 


-3.4422 


-3.9268 


-4.2875 


-4.7055 


-4.9258 


3s 










-0.1178 


-0.2758 




T=0.1 












Z 


3.0437 


3.0053 


3.0003 


3.0000 


does not 


does not 


2s 


-5.4903 


-6.1444 


-6.6282 


-6.9811 


converge 


converge 


2p 


-2.7462 


-3.4008 


-3.8853 


-4.2393 


due to 


due to 


3s 










3s level 


3s level 



effect and the onset of the decrease of R. Some authors 
have interpretted the plateau in the resistivity as indi- 
cating a situation where the mean-free path has become 
as small as possible ( as in the Joffc-Regel rule). We 
disagree with this explanation of the existence of a re- 
sistivity plateau in these cases in terms of a saturation 
of the mean-free path. Quantitative agreement is pro- 



vided with a very different picture of what happens in 
the plasma. As the material is heated, its resistivity in- 
creases just as in a metal, due to the increased avail- 
ability of a strip of states (of width ksT) at the Fermi 
surface for the "phonon-like" scattering to take place. 
However, as kgT increases, the number of current carry- 
ing electrons also increases, compensating the resistivity 
increase. During this process the chemical potential n of 
the plasma electrons begins to decrease, and a temper- 
ature T M= o is reached when /i passes through zero and 
towards negative values. Then the Fermi sphere is gets 
broken down and from then on, all the electrons, and 
not just those near the Fermi surface, begin to conduct. 
The plasma is essentially classical. The resistivity de- 
creases as the temperature increases. We are in fact in 
the Spitzer-like regime. The plateau defines the transi- 
tion to the Spitzer-like regime. The Fermi energy of the k 
= 0.1 case is small and it is already behaving like a clas- 
sical plasma. In reality the picture is more complicated 
since the Fermi surface changes not only because of the 
temperature, but also because the ionization of bound 
electrons changes the value of Z. This pushes the value 
of T M= o, an d the onset of the plateau to higher temper- 
atures than in a model with constant Z. It is easy to 
allow for this in a calculation of T M= o and confirm that 
the resistivity R(T, k) given in our table (and in Ref. 
is consistent with this picture. 



E. Contribution to the resistivity from 
electron-electron scattering. 

Discussions of the electrical resistivity of plasmas 
sometimes contain allusions to the e-e contribution to the 
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electrical resistivity. However, the electron-current oper- 
ator j — (e/m)k commutes with the electron-electron 
interaction Hamiltonian H PP . 



model where a model dielectric function is used. Thus, 



jH ee - H ee j = 



(17) 



This shows that the current is conserved under the e- 
e interaction. Hence electron-electron interactions can- 
not contribute to the resistivity arising from the elec- 
tron current. However, the e-e interaction has an indi- 
rect effect since it screens the electron-ion pseudopoten- 
tial V e -i(q). That is, electron-ion vertices must occur 
in all diagrams which contribute to the resistivity. In 
perturbation-theory approaches to the conductivity or re- 
sistivity, it is quite easy to get a contribution to R from 
e-e scattering alone, if the theory is carried out only to, 
say, second order. This means, if an all-oder calculation 
were done, the higher order corrections would exactly 
cancel the low-oder result for the e-e scattering contri- 
bution. Electron-electron interactions contribute to the 
resistivity of solids where the periodic potential generates 
Umkapp scattering. But this is not the case in plasmas if 
they can be considered uniform to within a length scale 
significantly larger than the mean free path. Classical 
transport calculations for systems with gradients (i.e, no 
translational invariance) is well known[26j. The quan- 
tum calculation is also well known in transport across 
heterostructures [22| . 



III. DYNAMIC CONDUCTIVITY. 

If we apply a field E cos(u;t) to the system, say using 
a light probe of frequency to, then the polarization of the 
medium is described by the polarization function H(q, ui) 
which is directly connected with the transverse dielectric 
function s(uj) and the dynamic conductivity cr(w). The 
wavevector q of the photon is nearly zero and may usually 
be omitted. The real part o~ r =tfto~{uj) at lo — > reduces 
to the static conductivity a that was already discussed. 



47TU> 2 



+ 1 



47T , , 

— n M 



(18) 
(19) 



If the effect of interactions and bound states is small, 
the conductivity of "free" electrons driven by the field 
Eoexp(itut), and damped by scatterng at ion centers is 
well approximated by the Drude model. It uses a relax- 
ation time To, (or a damping parameter 70) independent 
of the frequency. In partially ionized systems, or when 
there are interband effects in solids, it is necessary to in- 
clude bound-free and bound-bond contributions to a(to). 
A useful practical form is the extension of the Drude 



a r (u>) = 



■ l£i 



E 



UJEi(uj) 



-in 



(20) 
(21) 

(22) 



Here the pure free-electron Drude term uses the damp- 
ing parameter 70 and a plasma frequency oj p q, while 
other processes are modeled by a finite set of oscillators 
parametrized by 7^ and lu p \. This is a form useful for 
fitting experimental data since reflection and transmis- 
sion experiments could be used to extract best fit values 
of 7a,w p a, subject to the sum rule (f-sum rule): 



/ duja r (uj) = 
J ° \=0 



uj p \ = Tin tot e 2 /m e (23) 



Here ntot involves the free electrons (given by the average 
charge Z) as well as the electrons occupying the localized 
atomic states participating in the optical transitions. Un- 
like in liquid-metal studies, these are not readily available 
for warm dense partially degenerate plasmas. Further, 
the electron populations in the bound and free states are 
linked by the Fermi distribution associated with the elec- 
tron temperature and the chemical potential. Hence this 
LTE (local thermodynmaic equilibrium) type constraint 
also should be imposed in fitting the experimental data 
to the model dielectric function. Thus the most fruitful 
approach is to use the values of the mean ionization Z, 
occupation numbers, electron chemical potential fj, e , j\, 
e\ etc., provided from the density- functional NPA calcu- 
lation in setting up the fitting process. The simple Drude 
equation with a fixed value of 70 satisfies the f-sum rule. 
It can be shown that if we define n(w) in terms of the 
non-interacting polarizability H°(lj), by the form: 



n(w) = n°(w)/[i + G(w)n°(w)] 



(24) 



where G(u>) is called a "local-field" correction, then the 
sum rule is satisfied if it holds for n°(o>). Thus, if H(oj) 
can be calculated, then we can calculate the dynamic 
conductivity from it. 

The effect of corrections to the simple Drude term aris- 
ing from interband terms in solid Aluminum, were ex- 
amined many years ago by Dresselhaus, Harrison, and 
by Ashcroft and Sturm. 28]. They showed that there 
are band to band transitions near ^.05 and 0.15 of the 
Fermi energy Ep (which is ~ 12 eV. for normal density 
Al). The optical conductivity of normal-density liquid 
Al slightly above the melting point has been measured 
by Miller[2!| and shows a simple Drude form. A dis- 
cussion of liquid metal data including Cu, Ag and Au 
is found in Faberpjoj]. The deviations from the Drude 
form seen in solid- Al arise from the splitting of degen- 
erate bands due to the crystal potential, and to normal 
interband transitions. Benedict et al. have argued that 
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the thermal broadening of the electron self-energy is by 
itself sufficient to "wash out" these solid-state effects, 
even if the ion lattice remained intact |3lj. as may be the 
case in short-pulse laser generated plasmas. The effect of 
electron-electron as well as ion-electron interactions on 
the electron self-energy at finite-T was also discussed by 
us for the hydrogenic case [2(3 

When we consider warm dense ^M-plasmas, the Kohn- 
Sham eigenvalues can be used to assess if the driving 
field co could excite bond-bound or bound-free processes. 
As seen from Table [H] the 2s levels for the compressed 
systems (k =0.1, 0.5, 1, 2, 4) occur between -7.4 au to 
-5.5 au., while the 2p level ranges from -4.6 au. to - 
2.7 au. Hence, the Drude formula would be reasonable 
for these systems and for standard optical probes. The 
situation is quite different for the k =0.25 case. Here 
the 2p level ranges from ~ -4.5 au. to -5.5 au., but the 
3s level is very close to the ionization threshold. At T 

— 8 eV, the 3s level is at -0.165 au, and rises to -0.105 
au. near T= 3 eV, and then completely disappears for 
T below ~ 2.5 eV. Since a 300-400 nm optical probes 
corresponds to about 4-3 eV, it is clear that such probes 
would show deviations from the Drude conductivity for 
the k =0.25 case. When we go to low temperatures, the 
z = 3 ionization is very stable, while the 3s level creates 
the presence of z = 1,2 ionization states. In the case of n 

— 0.1, we have two shallow states , viz., 3s and 3p. The 
conductivity derived from the model dielectric function, 
i.e., Eq. 1201 would be relevant to the representation of 
experimental data for such systems. In the following we 
look at first-principles calculations. 

A. Perturbational approach to the dynamic 
conductivity. 

In a theory of a(co) we arc in effect looking for the po- 
larization function H(co). If the potential Vi e (q) is weak, 
(no bound states) diagrammatic methods can be used. 
However, the second-order expression in the screened in- 
teraction is essentially the only one that is tractable for 
realistic potentials. A number of such quasi-second order 
results exist in the literature, derived using various meth- 
ods. An old result, due to Hopficld, holds if the structure 
factor is essentially unity |32^. 

= "tAt / 7^3^l^ e (g)| 2 3[-^] (25) 
{Airco) z J (2nY e(q,uo) 

A more complete result, including the contribution from 
the dynamic structure factor of the ions can be written 
down, using an approach similar to that given by Mahan 
for phonons. Ropke et al. have also discussed second- 
order expressions within the Zubarev approach which is 
most suitable for short-ranged potentials free of bound 
states. [34j The approach used by Mahan is more inter- 
esting and may be used for long-range potentials. In fact, 
even in the two-temperature case where the ions are at a 
temperature Ti, while the electrons are at a temperature 



T e , it is easy to show using the Keldysh technique (as- 
suming that the probe frequency to does not overlap with 
core transitions) that the frequency dependent collisional 
relaxation time r{co) is|2E|. 

r (co) = -(u p u) J ~^1 Z J ^ (26) 

^[Xu(q, ^MXeefo v + w)]AJV(T e , Ti, V, CO) 

where we have set 

AN{T e ,T it v,u) = [N{P e {y + co)} - N{^(u)}}. (27) 

Here N{j3i(v)} is a Bose factor giving the occupation 
number of density-fluctuation modes at the temperature 
Ti and energy v. For the equilibrium situation we simply 
set Ti — T e . Then, as co — > 0, this equation can be 
shown to reduce to the inverse collision time used in the 
Ziman formula. In Eq. 1261 the electron response Xee and 
the ion-response \u mediate the energy and momentum 
exchange between the two subsystems. The imaginary 
parts of the response functions are related to the dynamic 
structure factors by a relation of the form: 

S(k,co) = -±coth[^}Z X (k,co) (28) 

In the case of the electrons, x(^> w ) can be written down 
in terms of the Lindhard function x (k, co) and the local- 
field correction G(k,co). In the limit where the elec- 
trons becomes classical, the Lindhard function simply 
reduces to the Vaslov function. The dynamic struc- 
ture factor of the ions can also be modeled in a similar 
fashion. [3^ Clearly, if the probe frequency co is smaller 
than the electron plasma frequency, then we may make 
a static approximation for G(/c, co). Since normal-density 
Al (and even Au, depending on the compression) have 
high plasma frequencies, this is a good approximation 
for dense Al or Au. However, the static G(k, co = 0) 
used should be such that at least the compressibility sum 
rule is satisfied. Also, we can introduce an co dependent 
model G(k,co) with G(k, co)= G(k,co — 0) for co < cop, 
and G{k, co) fitted to the high-frequency moment sum 
rules for co > cop . In practice it is not known how to sat- 
isfy all the sum rules. The static pair-distribution func- 
tion g{r) recovered from the imaginary part of such a re- 
sponse function should also agree with the known g(r) of 
the electrons at that density. The successful calculation 
of the g(r), S(k), and the G(k) in a consistent manner 
for electrons at strong-coupling and arbitrary tempera- 
tures (i.e, arbitrary degeneracies) was presented recently 
in Ref. |3| 

Although an extension of the above equation to take 
account of bound-free and bound-bound electron pro- 
cesses can be written down "by hand" , it is not easy 
to provide a systematic development. More fruit- 
ful approaches are to use density-functional molecular- 
dynamics simulations or time-dependent density func- 
tional theory. We consider these below. 
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B. Conductivity via the Kubo-Greenwood formula 
and molecular dynamics. 

Another approach to the conductivity a(u>) is to use 
molecular dynamics to develop the ionic-liquid structure, 
while retaining a DFT approach only for the electronic 
structure HUHlli^. The ion subsystem is modeled with, 
say, typically 32-256 atoms in a simulation box of volume 
Vb which is periodically repeated. An externally con- 
structed pseudopotential is used and an ionization model 
is assumed. The required number of electrons based on 
the ionization model is placed in the box. The ions are 
held at some fixed ionic configuration {Ri} and the Kohn- 
Sham electronic wavefunctions ipf(r, {Ri}) and energies 
e/ are computed. Given the size of the system, it is not 
practical to do more than a few k-points; usually only 
one fc-point, e.g., the T point is computed. The conduc- 
tivity for the given configuration, and for the selected 
k point with weight W(k, {Ri}) is estimated using the 
second-order Fermi golden rule formula. 

a(cj,{Ri},k) = -—J2\ W ( R ^ k )<^f\p\^9>\ 2 
fa 

(n f (k) - n g (k))S(e f (k) - e g (k) - u) 

This is in fact the Kubo-Greenwood (KG) procedure 
(sampled with one k-point) that one may use for a crys- 
talline solid. A Kohn-Sham form is used instead of the 
many-electron eigenstates. The occupation factors U/(fc) 
are also the one-particle occupations at the appropriate 
temperature. The energies used are the LDA eigenval- 
ues without corrections. Even a crystalline structure has 
its phonon modes and their effect is ignored in the en- 
ergy denominators. Similarly, the self-energy effects from 
the electron-electron processes are also ignored, although 
they can be quite large[2(| EJ. Only the umklapp pro- 
cesses associated with the reciprocal vectors of the sim- 
ulation box contribute to the static conductivity. This is 
because the wavefunctions are eigenstates of the "crys- 
talline" structure {Ri}- The actual conductivity has to 
be obtained by taking a configuration average which re- 
quires the Hclhmoltz free energies for all crystal configu- 
rations {Ri}. In effect, the method uses simple LDA- 
DFT for the electrons, but abstains from using DFT 
for the ions and carries out a detailed MD evaluation 
of the liquid structure. We believe that this is unnec- 
essary, especially for systems where the spherical sym- 
metry of the plasma is a statistically reasonable assump- 
tion. The work of Kown et al. 8] on strongly-coupled H- 
plasmas using these methods, and their comparison with 
our work is an example of this. In Fig. [3] we compare 
the ion-ion pair-distribution functions obtained from our 
NPA+ HNC+bridge type procedures with available sim- 
ulations from Silvestrelli et al 37j, and from Levesque et 
al|39|. A simulation box of 125 atoms implies only 5 
atoms per dimension and hence even the central atom 
feels only two atomic shells around it. The calculation of 
the u> — ► limit needed to obtain a static conductivity 
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FIG. 3: Ion-ion pair-distribution functions of Al from this 
work compared with those obtained by Silvestri|3'7^ . The in- 
set shows the comparison at normal density and at the melt- 
ing point. Here no data from ref. (37J is availble, but we use 
accurate molecular-dynamics simulations (Levesque et al[3f|) 
for comparison. 



is also quite difficult, and one approach is to assume the 
validity of the Drude form and fit a free-electron Drudc 
form to achieve this. The Kubo-Greenwod type MD pro- 
cedure would nevertheless provide useful complementary 
results for comparison with our two-component DFT ap- 
proaches, and would be of much interest in studying low- 
temperature systems with a tendency to covalent bond- 
ing and clustering. However, these methods bring in a 
number of approximations of their own and need cau- 
tious reconsideration. In this context, an interesting test 
would be to examine expanded liquid-Hg^3| using the 
Kubo-Greenwood MD approach. 



C. Dynamic conductivity from time-dependent 
density functional theory. 

Time dependent DFT provides a convenient approach 
to the calculation of the dynamic conductivity of in- 
teracting systems. The TDFT formulation relevant to 
dense plasmas [lOj includes dynamic screening and cou- 
pling to ion-dynamics in a computationally convenient, 
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eigenstates of the plasma. Then we have 

U(?,u) = U ext (f,u;) + V^ d + V^ d 



i 11(1 



v xc , = 



5n{r, w) 



\r — r 
dV xc (r : uj) 
dn(f, w) 

dr'xi?-, r'\u>)U(r', u>) 



5n(f, to) 



(30) 
(31) 

(32) 

(33) 



Here V^ d (r) is calculated from the gradient dV xc /dn 
evaluated at the density no(r) in the unperturbed neu- 
tral pseudo atom. Since the form of the time-dependent 
exchange-correlation potential V xc (r,u>) is still not es- 
tablished, most implementations use the static V xc [n(r)) 
of ordinary density functional theory. The above set of 
equations have to be solved self consistently to obtain the 
total perturbing potential U(r,oj). In effect, the dipole 
operator, or equivalently, the momentum operator of the 
scattering electron is replaced by a space and time depen- 
dent quantity which enters into the polarizability. Given 
the spherical symmetry of the system, the total perturb- 
ing potential has the form 



U(?,u) = -(l/2)EU u (r)Y°(r/r) 



(34) 



FIG. 4: Dynamic conductivity (au.) of Al at a compression 
of 0.25 at T =2.5 eV and 10 eV. The Drude conductivity 
is modified by the presence of the 3s level. Time-dependent 
DFT calculation greatly weakens the effect. 

self-consistent manner. The main consequence of the 
TDFT formulation is to replace, say, the dipole matix 
elements < i\r\j > between states i,j by a dynamic form 
< i|F(ci;)|j > where the modification of the driving field 
by the response of the system is taken into account in 
a self-consistent manner. Consider a weak external field 
E ex t(t) — Ecos(u>t). This corresponds to an external po- 
tential: 

U ext (r,t)=er-E ext (t) (29) 

The dipole form of the interaction is used since one of the 
objectives is to include the corrections arising from the 
presence of bound states. However, the dipole form of the 
matrix elements can be easily replaced by the momentum 
or acceleration formulation when needed. We assume 
that the electric field is directed along the z-direction, 
and suppress vector notation for the field unless needed. 
The external potential induces an electron density fluc- 
tuation 6n(r) which in turn generates corrections to the 
Coulomb and exchange-correlation potentials. Since the 
linear absorption coefficient (or optical conductivity) is 
the object of our study, 8n[r) etc., can be written in terms 
of the electron response function % (r, r'|o;) which can 
be approxmiately constructed|l0j from the Kohn-Sham 



Here Y® is a spherical harmonic. If we ignored these in- 
duced fields, the conductivity of the system can be writ- 
ten as: 

Sfla (u) = —3 V * Ni " 1 <HrV>| 2 (/. -Uu>) 
anh ^ 3 oj + e u — e' + iS 

(35) 

Here all the quantities on the right of the summation 
are in atomic units. The e 2 /(aoh) factor is the atomic 
unit of conductivity. Ni is the number of ions per unit 
volume, and we neglect the effect of ion-ion correlations 
in the bound-free and bound-bound processes. (It can 
be shown that these contribute mainly to the width of 
the transition by broadening the levels) Also when v, v' 
— k,l,m and k',l'm' for free-free transitions, the dipole 
matrix element is replaced by the momentum form, i.c, 
| < v\r\v> > | 2 = | < v\V\v' > \ 2 /u 2 . 

The conductivity expression cr°(w) is known to be par- 
ticularly inadequate when treating b-f and b-b processes, 
unless the initial bound state is a deep lying (e.g, Is) 
state^2j- Hence we need to include the induced fields 
in treating under-dense plasmas where there is a bound 
state close to ionization. This involves the use of the 
total perturbing potential, U(f,u>), rather than the ex- 
ternal potential r ■ E in constructing the matrix elements 
contained in the conductivity expression: 

I < v\U u {r)Y?(?/r)\v' > | 2 * (/„ - f nn ,)6(w + e v - k 2 /2) 
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The effect of level broadening can be included in the 
above expression by replacing the delta function S(lo + 
e„ — k 2 /2) by a form containing the self-energy correc- 
tions to the single-particle levels, as in Grimaldi et al, 
where self-energies were calculated. However, for our 
present purpose, we use the r used in the Drude cal- 
culation to provide a broadening parameter. In fact, we 
have shown in Ref. 20] that the self-energy contribution 
of ion-electron scattering to level broadening is identi- 
cal to that given by the Ziman formula. Thus, setting 
7 = 1/t, we replace the 5-function in the above equation 
with a Lorentzian. 



S(lu + e v - k 2 /2) 



7/tt 



(w + e v - fc 2 /2) 2 + 7 2 



D. Some numerical results. 

In this section we present results for the dynamic con- 
ductivity of some Al-plasmas to illustrate the effect of 
the scattering events where shallow bound states mod- 
ify the Drude-like conductivity. These "bound states" 
are really Kohn-Sham eigenvalues and hence their val- 
ues may need improvement by evaluating the corrections 
using a Dyson equation. [2(j The first step of the calcula- 
tion is the solution of the NPA model to obtain the Kohn- 
Sham basis set {4>{t)u} for the given electron density and 
temperature. The codes necessary for the NPA calcula- 
tion, the resistivity calculation, as well as the plasma 
conditions (EOS) that go into the resistivity or conduc- 
tivity calculations may be accessed via the internet at our 
website. 18]. Some results for the shallow bound states 
present in AZ-plasmas at T=2.5 eV and 10 eV, for com- 
pressions k = 0.25, and 0.1 are given in table II. In the 
T=10, k = 0.1 case we have two shallow bound states, 
viz., 3s and 3p. In Tabic ITTT1 their occupation numbers 
and the interacting chemical potential needed for calcu- 
lating the Fermi factors are given. In calculating the ma- 
trix element | < nlm\r\kl'm' > | 2 , the correct density of 
states H(e) for the continuum state \kl'm! >, normalized 
in a sphere of radius R, with R — > oo should be included. 
That is, 
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FIG. 5: Dynamic conductivity (au.) of Al at a compression 
of 0.1 at T =2.5 eV and 10 eV. The Drude conductivity is 
modified by the presence of shallow 3s at T = 2.5 eV, and 
also the 3p level at 10 eV. The effect is reduced in the TDFT 
calculation. 



be regarded as highly provisional results only; in fact, 
the treatment of shallow boundstates using DFT is itself 
questionable since the Kohn-Sham eigenvalues are known 
to be a poor approximation to the actual excitation spec- 
trum. 



k n R-nl/2 



dndk 
H(£) = dkd~e 

The calculations are very simple if (a) we ignore the 
Zangwill-Soven type TDFT effects arising from the re- 
action of the system which act to modify the external 
field, (b) if we ignore the phase shifts Ski, and replace 
the boundstates by hydrogen-like states with the correct 
Z = Z. The final results depend on the level broadening 
parameter 7 used in the calculation. 

In fig0]and Fig[5]we show the conductivity arising from 
the presence of the 3s and 3p levels, as well as the Drude 
term, and the modification when TDFT is included in an 
approximate manner. These TDFT calculations should 



bound-bound processes. 



Clearly, the presence of a partially occupied 3s and a 
3p state with an energy separation of the order of 0.46 au. 
would lead to a contribution from b-b processes at around 
12 eV. This contribution is easily included in the calcula- 
tion under the simplifying assumptions that we noted be- 
fore. However, this is a relatively sharp "line" resonance 
which is expected to undergoes significant modification 
when the time-dependent effects are taken in to account. 
We have not included it in our figures. 
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TABLE III: details regarding the n=3 shallow states in un- 
derdense Al at T=10 eV and 2.5 eV. Full occupation is when 
Occ is unity 



K 


0.25 


0.25 


0.1 


0.1 


T 


2.5 eV 


10 eV 


2.5 eV 


10 eV 




0.0393 


-0.5538 


-0.6652 


-0.9589 


Occ(3s) 


0.7444 


,0.2217 


0.6850 


0.1455 


Occ(3p) 








0.0843 



IV. CONCLUSION. 

The first-principles calculation of the dynamic conduc- 
tivity of warm dense matter may be conveniently car- 



ried out within the framework of multi-component den- 
sity functional theory. The static calculation (NPA etc) 
provides the Kohn-Sham basis set, phase shifts, pseu- 
dopotentials for constructing ion-ion pair potentials and 
structure factors. These immediately provide results for 
the static conductivity. Further, the energy-levels and 
occupation numbers obtained from the Kohn-Sham NPA 
solution can be the starting approximation for a time- 
dependent density functional calculation of the optical 
conductivity. This proceeds in much the same way as the 
optical absorption cross section. Details of such a time- 
dependent calculation, based on the method of Zang- 
will and Soven, may be found in Grimaldi, Lecourte and 
Dharma-wardana. lOi 
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